Fluctuations in a fluid with a thermal gradient 
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We report measurements of the temperature fluctuations in a horizontal layer of a pure fluid 
confined between parallel plates and heated from below. Consistent with earlier work we found that 
the structure factor Sriq) (the square of the modulus of the Fourier transform of the temperature 
field) was consistent with Sriq) ~ q~'^ for large horizontal wavenumbers q. As q decreased, S{q) 
increased less rapidly, passed through a maximum, and then approached zero. The results agree 
qualitatively, but not quantitatively, with the predictions of Ortiz de Zarate and Sengers. 

PACS numbers: 05.40.-a, 47.20.-k 



The enhancement of fluctuations in a fluid due to non- 
equilibrium conditions has long been an important topic 
in statistical mechanics. Even for the more specific case 
of a fluid with a thermal gradient this problem has a 
considerable history. [1-3] For large wavenumbers q per- 
pendicular to the gradient it was shown two decades ago 
that the structure factor S{q) (proportional to the square 
of the modulus of the Fourier transform of density fluc- 
tuations) should vary as [4-8] This was confirmed 
experimentally by small-angle light scattering. [9-12] It 
was soon realized theoretically that the small-g diver- 
gence saturates at a constant value as q vanishes because 
of gravity. [13,14] This saturation effect was indeed ob- 
served by light scattering at very small angles from binary 
mixtures where concentration fluctuations arc expected 
to play a role similar to temperature fluctuations in a 
pure fluid and give a much stronger signal. [15,16] For 
a sample confined between parallel plates an additional 
saturation effect is expected to come into play because 
fluctuations orthogonal to the plates are cut off by the 
finite plate spacing d. [17,18,2,3]. With decreasing q this 
I latter effect is expected to lead to a maximum of S{q) , 
with a subsequent small-q regime where S{q) ^ q^ . How- 
\ ever, to our knowledge this maximum and the decrease 
at even smaller q was not accessible to prior experiments. 
' Very recently Ortiz de Zarate and Sengers (OZS) [3] cal- 
culated quantitatively the shape of S{q) over the entire 
I range of q for realistic no-slip top and bottom boundary 
conditions, using the linearized Boussinesq equations of 
' fiuid mechanics [19] and a low-order Galerkin approxi- 
mation for the fluctuations in the vertical direction. 

I Here we report measurements of Sriq) of temperature 
fluctuations which are related to S{q) by a thermody- 
namic factor. We covered the range 1 ^ q ^ 6.5 where 
q — 2tt/X with A the horizontal wavelength of the fluc- 
tuations in units of d. We used a shadowgraph method 
and a thin horizontal layer of fluid heated from below. 
The sample was compressed SFg on the critical isochore 
above the critical temperature Tc = 45.567°C where fluc- 
tuation effects become exceptionally large, and where the 
shadowgraph method becomes very sensitive because the 



refractive index n has a large temperature dependence 
{dn/dT)p. We found that Sriq) had the expected max- 
imum in the range 3 ^ g 4, and that it decreased 
for larger g in a way consistent with S't('z) ~ q~'^ as 
predicted and previously observed. Below the maximum 
Sriq) tended toward zero with decreasing g in a manner 
consistent with the q^ dependence predicted by OZS. Al- 
though our results agree qualitatively with the prediction 
[3] by OZS, the agreement is not quantitative. The rea- 
son for the difference between theory and experiment is 
not clear at present. On the one hand it might be found 
in the fact that the calculation is based on a low-order 
Galerkin approximation. On the other, it could be due to 
nonlinear interactions between the fluctuations and the 
use of linearized equations of motion in the theory. 

The apparatus was described in detail elsewhere. [20] 
The top and bottom plates of the cell each consisted of 
a 0.318 cm thick sapphire plate with a diameter of 10 
cm. [21] A silver layer of thickness 0.9 fim was evapo- 
rated on the top surface of the bottom plate to provide 
a mirror for the shadowgraph. An aluminum plate of 
equal size containing thermistors and a film heater was 
pressed from below against the bottom sapphire. The 
separation d between the sapphire plates was fixed by 
a porous paper ring with an inner (outer) diameter of 
2.5 (3.3) cm. It was measured interferometrically [20] 
to be d = 34.3/im. Inspection with an expanded laser 
beam revealed circular Newton rings with their centers 
close to the center of the paper ring. The actual area 
sampled by the shadowgraph was near the center of the 
Newton rings where convection rolls first appeared when 
AT was increased above the onset at ATc, and had a 
size of 1.3 X 1.3 mm. Assuming a parabolic profile of the 
bowing of the top/bottom plates, we estimate that d was 
uniform to ±0.017/im or 0.05 % over the sampled area. 
About half the temperature difference across the cell was 
across the SFg layer. The measured AT was corrected 
accordingly. This introduced a systematic uncertainty of 
AT of about 10%. Whereas this primarily leads to an 
overall factor of order unity in the comparison with the- 
ory, it makes a sufficiently accurate a priori adjustment 
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of the mean sample temperature T to the value T{pc) on 
the critical isochore difficult. Thus we determined T{pc) 
experimentally. Figure 1 illustrates the phase diagram 
near the critical point in the T — p plane. The pressure 
was held constant within 0.0005 bar at 38.325 bars, cor- 
responding to the solid line. In order to find T{pc), the 
fluctuation power Pj at constant AT was measured as a 
function of T by the shadowgraph method. The value of 
T at which Pj had a maximum was taken to correspond 
to T{pc). Later, when AT was changed in the experi- 
ment, T was held constant at T{pc). Thus the imposed 
AT caused the sample density to vary vertically between 
the hot and cold plates as illustrated by the heavy solid 
line in the figure. 
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1. Temperature-density plane near the critical point 
Dashed line: coexistence curve separating liquid and 
Vertical dotted line: critical isochore. Solid circle: 
critical point Tc = 45.567° C, Pc = 37.545 bar, pc = 0.742 
g/cm^. Solid line: isobar P = 38.325bars. Heavy solid 
line ending in the two circles: density range spanned for 
AT = 0.220°C. 



The shadowgraph apparatus was a modified version 
of one described before. [20] We used a digital camera 
with 12 bit resolution. [22] The light source was a light 
emitting diode [23] behind a nominally 15/im diameter 
pinhole. A coUimating lens with focal length / = 5 cm 
was used. The shadowgraph transfer function T{q,zi) 
was measured by determining the fluctuation intensity 
at constant g as a function of the optical distance zi 
[20,24]. We typically took N = 1000 uncorrected im- 
ages li of size 272 x 272 pixels at a given AT and used 
their average Jq as a reference to compute the signal 
7j = /j//o — 1. The square of the modulus of the Fourier 
transform Si (q, AT) of li was used to calculate the struc- 
ture factor S'T,i(q, AT) of the temperature field of each 
image. [20,24] The average St{ci, AT) of 5T,i(q, AT) over 
all TV images was used to compute the azimuthal integral 
ST{q,AT). 

For AT = the experimentally determined 5t((7, 0) 
was dominated by instrumental noise. Thus we derived 
the difference SSriq, ^T) = Sriq, AT) - Sriq, 0). It fol- 
lows from the work of OZS that this quantity is given 
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The factor 5^ gives the strength of the temperature fluc- 
tuations in the equilibrium system and should not depend 
on q or AT. It contains a number of thermodynamic pa- 
rameters and Boltzmann's constant. Its estimate from 
known thermodynamic data is not highly accurate, and 
we chose to treat it as an adjustable scale factor. Accord- 
ing to OZS the strength of the nonequilibrium enhance- 
ment is given by [25] 
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Here Cp is the heat capacity at constant pressure and 
T the absolute temperature. The parameter a = v/Dt, 
with V equal to the kinematic viscosity and Dt equal 
to the thermal diffusivity, is the Prandtl number and 
R = apgd^ AT / {uDt) , with ap equal to the isobaric 
thermal expansion coefficient and g the acceleration of 
gravity, is the Raylcigh number. Within the low-order 
Galerkin approximation used by OZS the dependence 
of the nonequilibrium enhancement on the horizontal 
wavenumber q is given by 
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In Fig. 2 we show experimental results for SSxiq, AT) 
for several AT. Note that all of the AT are well below 
the value for the onset of convection AT^ = 0.440°C. At 
large q the data show the rapid decrease with increas- 
ing q which corresponds to SSxiq, AT) ~ q~^ and which 
had been observed by light scattering. As q decreased, 
the q~^ dependence saturated because fluctuations in the 
vertical direction were inhibited by the boundaries. At 
even smaller q the structure factor had a maximum and, 
with decreasing q, vanished. This latter region had not 
been observed before in experiment, but was predicted by 
OZS. The maximum of SSriq, AT) occurred at a value 
of q which corresponds roughly to the inverse of twice the 
plate spacing. All the qualitative features are consistent 
with the theoretical prediction. 

At a quantitative level the agreement with theory is 
not perfect. We fitted the prediction Eq. 1 to the data, 
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adjusting S]^ separately for each AT. The results are 
given by the solid lines through the data sets in Fig. 2. 
Clearly the theoretical curves are much wider than the 
data suggest. On the other hand, the results for S'U were 
only mildly dependent on AT, suggesting that the de- 
pendence on AT contained in Eq. 2 is about right. 
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FIG. 2. Experimental results for the structure factor 
SSriq, AT) as a function of the wavenumber q for several AT. 
From top to bottom the results are for AT =0.189, 0.142, and 
0.095 °C. For comparison, the onset of convection occurs at 
ATc = 0.440° C. The lines are fits of Eqs. 1 to 6 to the data, 
adjusting only S%. 



curves, i.e. wc replaced R with an adjustable R in Eqs. 2 
and 5. [25] We show the resulting fits in Fig. 3. To better 
reveal the behavior at large and small q, the same results 
are shown on double logarithmic scales in Fig. 4. We note 
that the fits of these empirically modified theoretical pre- 
dictions are excellent for all three AT. Prom Fig. 4 we 
also see that the data agree well with the predicted 
and dependence at large and small q respectively. 
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FIG. 4. Experimental results for the structure factor 
SSriq, AT) a function of the wavenumber q as in Fig. 3, 
but on logarithmic scales. 
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FIG. 3. Experimental results for the structure factor 
SSriq, AT) as a function of the wavenumber q for several 
AT. The data are the same as those in Fig. 2. The lines are 
fits of Eqs. 1 to 6 to the data, adjusting the two empirical 
parameters fq and R discussed in the text in addition to S%. 



In order to obtain a better fit, we introduced two em- 
pirical AT-dependent adjustable parameters in addition 
to the overall scale factor S]^. First, we allowed an ad- 
justment of the length scale in Eqs. 5 and 6 by replac- 
ing q with fqq and least-squares adjusting fg. This has 
the potential to expand/contract the curves horizontally 
along the g-axis. Initially this was done to accomodate a 
possible experimental uncertainty in d which enters into 
computing q. However, this uncertainty allows only val- 
ues in the range = 1 ± 0.03 and is independent of AT. 
The second adjustment was made in the efi'ective value 
of R because this parameter influences the widths of the 
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FIG. 5. Experimental results for R/Rc (solid line) and for 
the fit parameter R/Rc (solid circles) as a function of AT. 
According to the results of OZS the two should be equal. 
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FIG. 6. Experimental results for the length-scale fit param- 
eter fq as a function of AT. According to the results of OZS 
fq should be equal to one. 
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FIG. 7. Experimental results for the coefficient as a 
function of AT. According to the prediction of OZS S% 
should be independent of AT. 

In Fig. 5 we show the results for R divided by the re- 
sult for Rc given by the Galcrkin calculation of OZS, 
i.e. by Rc = 1750. This parameter primarily influ- 
ences the width of SSriq)- Also shown, as a solid line, 
is R/Rc — AT/ATc based on the experimentally de- 
termined ATc = 0.440°C. According to the theory the 
two should be equal. As AT approaches ATc = 0.44°C, 
the fitted values of R/Rc approach the expected value 
R/Rc = 1, but at smaller ATc there is a large disagree- 
ment over a wide range. 

Figure 6 shows the length-scale parameter fq. It is less 
than the predicted value of unity over the entire range of 
the experiment. It can reasonably be extrapolated to 
fg = 0.97 at AT ^ ATc 0.440°C. Within experimental 
uncertainty this is consistent with fq(Rc) = 1- 

Finally, in Fig. 7 we present the results for S'^ which 
control the overall size of SSriq)- According to the the- 
ory this quantity should be independent of AT. The 
experimental results are indeed reasonably constant. 

In this Letter we presented experimental results for the 
non-equilibrium contribution 5St{(1, AT) to the struc- 
ture factor of a horizontal layer of fluid confined between 
parallel plates and heated from below. We found that the 
results have the qualitative features of the predictions of 
OZS. In particular, at large and small q we found the 
data to be consistent with 5St ^ q~^ and 5 St ^ q^ re- 
spectively, with a maximum of 6 St {q) at an intermediate 
value of q comparable to that given by the inverse plate 
spacing. To our knowledge the maximum, due to the con- 
finement, and the dependence at small q had not been 
observed before in experiments. A detailed comparison 
of the shape of SST{q) reveals quantitative differences be- 
tween experiment and theory. In order to obtain a fit to 
the data, we introduced two empirical parameters. One 
of them, an effective Rayleigh number i?, modifies the 
width of the curves. As shown in Fig. 5, R differs from R 
over the entire range of AT, with the difference largest 
at small AT. The other, fg, modifies the characteristic 
length scale. As seen in Fig. 6, it also differs from the 
prediction fq = 1 over the entire range, again with the 
largest difference at small AT. 

At this point we can only try to speculate about the 
reason for the difference between theory and experiment. 
Two possibilities come to mind. On the one hand, the 
low-order Galerkin approximation which was used may 
be inadequate. This could be tested by (numerical rather 
than analytic) higher order calculations. On the other, 
the neglect of nonlinear terms in the deterministic equa- 
tions of motion may not be appropriate for the conditions 
of the experiment. However, in that case we would have 
expected the disagreement between experiment and the- 
ory to diminish as AT decreased. As seen in Figs. 5 and 



6, just the opposite is the case. 
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